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Abstract 

The Falicov-Kimball model of spinless quantum electrons hopping on a 1-di- 
mensional lattice and of immobile classical ions occupying some lattice sites, with 
only intrasite coupling between those particles, have been studied at zero tempe- 
rature by means of well-controlled numerical procedures. For selected values of 
the unique coupling parameter U the restricted phase diagrams (based on all the 
periodic configurations of localized particles (ions) with period not greater than 16 
lattice constants, typically) have been constructed in the grand-canonical ensemble. 
Then these diagrams have been translated into the canonical ensemble. Compared 
to the diagrams obtained in other studies our ones contain more details, in particular 
they give better insight into the way the mixtures of periodic phases are formed. Our 
study has revealed several families of new characteristic phases like the generalized 
most homogeneous and the generalized crenel phases, a first example of a structural 
phase transition and a tendency to build up an additional symmetry - the hole- 
particle symmetry with respect to the ions (electrons) only, as U decreases. 
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1 The model and the method 



The work reported in this paper has been inspired by the recent results due to Freericks 
and Falicov |1|], Farkasovsky and Bat’ko 0, and Gruber et al. 0. It is also a continuation 
of our earlier work 0 . All the papers mentioned are concerned with various properties of 
the ground state phase diagram of the two-state spinless Falicov-Kimball model. This is a 
lattice model of itinerant quantum electrons and immobile classical ions that are coupled 
via a simple on site interaction. It is a simplified version of the model that was put 
forward in 1969 by Falicov and Kimball || to describe conductor-insulator transitions. 
In the one-dimensional case, considered in this paper, the Hamiltonian corresponding to 
the model has the form 

H a = J2 (~ a *x a x+i ~ <+i a x - Uw(x)a*a x ^ , (1) 

where A is a finite piece of the one-dimensional lattice with periodic boundary conditions 
imposed, a*, a x are the creation and annihilation operators of itinerant spinless electrons 
and w(x) = 0 or 1 is the occupation number of the ions at the lattice site x. The function 
A 9 x — > w(x), is called the ion configuration. Observe that the electrons can hop only 
between the nearest neighbour sites and the energy units have been chosen in such a way 
that the hopping rate is unity. Therefore in H\ there is only one independent parameter 
expressed in these units: the electron-ion coupling U. The total number of electrons 
N e = Yv a * x a x and the local occupation numbers w(x) of the ions (hence the total number 

x£A 

N t (w) = X] w(x) of the ions in the configuration w) are the conserved quantities. 

x£A 

The Falicov-Kimball model can be related to the celebrated Hubbard model ||, 0. 
It has however many interpretations of its own [[L]], j8|, interesting enough to motivate 

intensive investigations of this model. In one of the interpretations proposed by Khomskii 
||, Kennedy and Lieb p] and adopted here it is a model of crystalline formation. As 
proved by Kennedy and Lieb || and by Brandt and Schmidt || the effective interaction 
between the ions, induced by the direct interaction with mobile electrons, leads to the 
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formation of periodic structures of the ions (superlattices). Another interpretation that 
is sound from the point of view of solid state physics, is in terms of mobile s-electrons 
and localized /-electrons in rare-earth compounds. This makes the model relevant for 
understanding mixed- valence phenomena 0, j8|. 

For a given ion configuration w the Hamiltonian (|l[) is the second quantized form of 
the one-particle operator h(w ) whose matrix elements are 

[h(ic)]a;j/ txy Uw (x) S X y , (2) 

where 

+ = / 1 if y = x ± i ( o\ 

xy 1 0 otherwise. 

In the sequel by the spectrum of the configuration w we mean the spectrum of h{w). 

The states of the whole system are characterized by their ion configurations and the 
electrons’ wave functions. If the ion configuration w and the total number of electrons 
N e are specified, then in the ground state the electrons’ wave function corresponds to the 
so called Fermi sea. The energy of this state, E u (w, N e ), called the Fermi sea energy, is 
equal simply to the sum of all the N e lowest eigenvalues of h(w). Such states of the system 
are referred to as phases. Suppose that we are given the total number N e of the electrons 
and the total number of the ions. In the ground state subject to these conditions the 
ions are arranged in a specific configuration w, the ground state configuration, and the 
electrons form the corresponding Fermi sea. This state of the whole system is referred 
to as the (N^Nf)- phase. Its energy, E(N e ,Ni), called the ground state energy, is equal 
to E u (w, N e ). The largest eigenvalue of h[w) that contributes to E L (w,N e ), referred to 
as the Fermi level, is separated from the higher and lower parts of the spectrum of w by 
the gaps A u (w ) and A i(w), respectively. A phase that is the (N e , A*)-phase for some pair 
(N e , Nf) is said to appear in the canonical phase diagram or is referred to as stable phase. 

Our results refer exclusively to the infinite system. The thermodynamic limit is 
taken in such a way that the electron number per site N e / |A|, the ion number per site 
Aj/|A|, the Fermi sea energy per site E u (w, N e ) / |A| and the ground state energy per 
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site E(N e , Ni)/\A\ converge to the electron density p e , the ion density the Fermi sea 
energy density e u (w, p e ) and the ground state energy density e(p e ,pi), respectively. Sup- 
pose that the ion configuration of the (p e , -phase is w. If in the thermodynamic limit 
A u (w) — A i(w) > 0, then the (p e , Pi )- phase is called the insulator , otherwise it is termed 
the conductor. 

The possibility of studying directly infinite systems stems from the fact that for the 
periodic ion configuration w the Green function of the tight binding Schrodinger equation 
and consequently the density of states n(w), p e and e u (w; p e ) can be determined exactly. 
However the corresponding formulae depend on the band edges of the spectrum of w, 
which are zeroes of the polynomial whose order coincides with the period of w 0- First 
these zeros have to be determined numerically and after that the quantities of interest, 
like the energy density e u (w; p e ), can be obtained by a numerical integration. 

The hole-particle unitary transformations applied to the Hamiltonian (|I|) reveal the 
corresponding symmetries of the Fermi sea energy density e u {w\p e ). Thus, performing 
the hole-particle transformation with respect to the ions we obtain 



e {u) (w*i pe) = e { U) (tu; p e ) - U p e , 



- J-u) 



( 4 ) 



where w* is the adjoint ion configuration of w, w*(x) = 1 — w(x), while with respect to 
the electrons we get 

e (t7) (w; p e ) = e ( ~ u) (w; 1 - p e ) - Up^w). (5) 

Therefore we can reduce the range of p e values considered to the interval [0, 1/2] and we 
can also fix arbitrarily the sign of U. We choose U > 0, i. e. according to the ions 
and electrons attract each other at the same lattice site. Thus it is natural to refer to the 
system with p e = Pi as neutral. 

In this paper we study the properties of the Falicov-Kimball model at zero temperature 
by means of an approximate method, the method of restricted phase diagrams. It is the 
set of the states available for the system considered that is approximated (by a restricted 
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set of states), not the interactions. We concentrate mainly on the ion subsystem, that is, 
we try to determine how the ions are arranged in the ground state. This sort of approach 
has already been used, first in [|1T] , then in [j|l|], [|2|, J3j], Q, 



The method of restricted phase diagrams is straightforward in principle. In the canoni- 
cal ensemble, we are interested in here, some specified values are assigned to the densities 
p e and pi and the task amounts to minimizing the set of functions {e(w,p e ) : w € W Pi } 
over the set >V Pi of admissible configurations that satisfy the condition pfiw) = p*. The 
minimum function constitutes the required approximation to the ground state energy 
density e(p e , pfi . According to the thermodynamic stability conditions e(p e , pfi should be a 
convex function of (p e , pfi). Since the functions p e — > e{w, p e ) are convex and the minimum 
of a set of convex functions is not necessarily convex, to guarantee the thermodynamic 
stability it is necessary to construct the convex envelope of the minimum function. It 
is clear that the choice of W Pi is decisive for the quality of the approximation. We can 
choose W Pi to be for instance the set W PitP of all periodic ion configurations whose period 
does not exceed some p and pfiw) = pi. Then those points of the graph of the minimum 
function that coincide with the graph of the convex envelope can be associated with 
the periodic phases. The remaining points can be obtained as convex combinations of 
some other points belonging to the graph of the minimum function and therefore they 
are associated with mixtures of periodic phases. The canonical restricted phase diagrams 
of this kind have been obtained in |I|] and |2|. It turns out however that the resulting 
canonical restricted phase diagrams are highly unstable. There are two main reasons for 
lack of stability in these diagrams. First, they do not allow for mixtures of periodic phases 
with different ion densities. In many instances the periodic phase corresponding to the 
point (p e , p^ of the restricted phase diagram can be shown to have higher energy than a 
suitable mixture of two or three periodic phases whose ion densities are different from pi 
(cf. 0 and section 3). Second, the sets of admissible ion configurations used there are 
too restricted. 
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To overcome the stability problems mentioned and to avoid dealing with mixtures we 
have studied the system in the grand canonical ensemble and then translated the results 
into the canonical ensemble. In the grand canonical ensemble one is given the electron 
and ion chemical potentials (He, Hi) and the task is to find the minimum function of the 
set {/ (in; He, Hi) : w E W} of the free energy densities / ( w ; p e , Hi), where 

/ O; He, Hi) = e ([/) (in; He) ~ He Pe O; He) ~ Hi Pi( w ) (6) 

and W is a set of admissible configurations. The thermodynamic stability requires the 
minimum function to be a concave function of (p e , Hi)- Since the free energy / ( w ; He, Hi) is 
a concave function of (/i e , Hi) and the minimum function of any set of concave functions is 
concave, the thermodynamic stability condition is satisfied automatically. We have taken 
for the class W all periodic ion configurations whose periods do not exceed 16 lattice 
constants, typically, however we have also constructed some diagrams where the maximal 
period has been increased to 40. 

The canonical phase diagrams have a [/-independent inversion symmetry about the 
point (1/2, 1/2). It follows from ([Ij) and (j5|) that, if w is the ground state configuration 
of the ions at the point ( p e , Pi ), then the adjoint configuration w is the ground state 
configuration at (1 — p e , 1 — p*). 

2 Some remarkable phases and their ions’ configura- 
tions 

To each point (p e , Pi) the canonical phase diagram assigns a phase (more precisely a class 
of equivalent phases that consists of the phases related by the symmetries of the lattice), 
which is called here (p e , /?j)-phase or the stable phase. Since we are concerned mainly 
with the ion subsystem, we characterize the phases by their ions’ arrangements w and 
the electron density p e . From a few thousand of the periodic ion configurations that have 
been tested, each time a restricted phase diagram was constructed, only a few per cent, 
at most, happen to be the ions’ arrangements of the phases that appear in the phase 
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diagram. Many phases turn out to be the members of some infinite families of phases 
whose periodic ion configurations comply with peculiar rules and the ion and electron 
rational densities satisfy linear relations. 

There are of course only two translationally invariant ion configurations, the full con- 
figuration, where all the sites of the lattice are occupied (w(x) = 1) and the empty 
configuration with no ions at all (w(x) = 0). The phases with translationally invariant 
ion configurations are called the full phases and the empty phases, respectively. The 
unique band of the empty configuration extends from —2 to 2, while the band of the full 
configuration is obtained by translating the latter by —U. In the full or empty phases 
the system is a conductor, unless p e = 0 or 1. The spectrum of any ion configuration is 
confined to the interval given by the lower band edge of the full configuration, —2 — U, 
and the upper band edge of the empty configuration, 2. For U > 4 the two bands do not 
overlap and the gap separating them opens in the spectrum of any other configuration. 

From now on when we speak of periodic phases or periodic ion configurations we mean 
that their period is at least 2. 

The infinite families of the ion n-molecule most homogeneous phases [p/q] n are char- 
acterized by the so called ion n-molecule most homogeneous configurations of ions and 
the relation = np e , where n is an integer, p e = p/q and 0 < p e < 1 fn. In the context 
of the Falicov-Kimball model the ion n-molecule most homogeneous configurations have 
originally been defined in |l2j (the case n = 1) and in [[| (n > 2). Let p and q, p < q, be 
relatively prime integers. The unit cell of the ion n-molecule most homogeneous configu- 
ration, whose ion density is pi = np/q, has the length q and the arrangement of the ions 
is given by the solutions kj to the equations 

p kj = j mod q, j — 0, 1, . . . , np — 1. (7) 

The ions form clusters, the ion n -molecules, that consist of n consecutive sites occupied 
by the ions. The relation between the ion and electron densities means that in these 
phases there is one electron per ion n-molecule. The distances between any two nearest 
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neighbour clusters are either d or d + 1. Similarly the distances between next nearest 
neighbour n-molecules are either d! or d! + 1, and so on. In |3[] the ion 1-molecule most ho- 
mogeneous configurations have been called the atomic most homogeneous configurations. 
The following example illustrates the definitions introduced: for p = 2 and q = 7 the unit 
cell of [2/7] i phase is [* 000 * 00 ], the unit cell of [2/7]2 phase is [• • o o • • o], and the 
unit cell of [ 2 / 7)3 phase is [•••••• o]. 

The infinite families of the hole n-molecule most homogeneous phases [p/ q\* n are char- 
acterized by the hole n-molecule most homogeneous configurations of the ions and the 
relation 1 — Pi = np e , where n is an integer, p e = p/q, 0 < p e < 1/n. The hole n-molecule 
most homogeneous configuration is defined in a similar way as the ion n-molecule most 
homogeneous configuration, except that the ions are replaced by the holes. In this sense 
the star can be interpreted as the operation of taking adjoint. Clearly the adjoint of the 
adjoint phase is the original phase. In these phases there is also one electron per molecule 
but this time it is the hole n-molecule. 

In each family of the phases described so far the corresponding periodic ion arrange- 
ments can be viewed as composed of the same sort of molecules and there is one electron 
per molecule. There appear however in the phase diagrams other families of phases, where 
the corresponding periodic ion arrangements can be viewed as composed of the hole n- 
molecules, with n > 2, and of the composite molecules that consist of 1 ion, 1 hole and 
1 electron. In these phases there is also one electron per each composite molecule but 
the number of electrons per hole molecule depends on its size. Specifically for a pair of 
irreducible rationals p/q and h/k , such that 1/3 < h/k < p/q < 1/2 and k — 2h — 1, the 
\p/q\h/k phase has period q, the electron density p e = p/q, the ion density pi = kp e — h, 
and the positions of the ions in its unit cell are given by the solutions kj to the equations 

p kj = j mod q, j — 0, 1 , . . . , kp — hq — 1 . (8) 

The rationals h/k of \p/q\h/k phases that are stable in our restricted diagrams belong to 
a special sequence, which will be specified later. The [p/q]h/k phase can be viewed as 



a periodic arrangement of the hole fc-molecules with h electrons per hole molecule and 
composite molecules built of 1 ion, 1 hole and 1 electron. This definition of \p/q\h/k 
families of phases does not apply to the limiting case p/q = h/k = 1/2 whose definition 
is postponed to the text that follows. 

In the phase diagrams there appear also the families of phases that are adjoint to 
the [p/q\h/k phases. Every \p/q\h/k phase has its adjoint phase \p/q\* h / k , where the hole 
fc-molecules are replaced by the ion fc-molecules. The ion and electron densities of the 
\p/q]1/h phase satisfy the linear relation p^ = —kp e + k — h. 

One of the remarkable features of the families of periodic phases described so far is 
that their density pairs (p e , Pi) and their ions’ configurations can be obtained recursively, 
following a recursive construction of the so called Farey sequence. Therefore we can 
speak of the Farey sequences of phases like the Farey sequence of \p/q\ n phases or the 
Farey sequence of \p/q}h/k phases, etc. 

The Farey sequence F n of order n is the ascending sequence of irreducible rational 
fractions between 0 and 1 whose denominators do not exceed n, i.e. the fraction p/q , with 



p and q relatively prime integers, belongs to T n ifO<p<g<n |H|. Having any two 
consecutive elements p/q and p'/q 1 of jF n , the ’’parents”, we can form their ’’descendant” 

p"/q ", 

P + P ' ( 9 ) 



P 



q + q r 



which satisfies the inequalities 



p p p 

q q" q' 



( 10 ) 



This suggests how to construct T n recursively. To build up T n we start with which is 
simply the pair of the ’’first parents” {0/1, 1/1}, and add to it the descendant 1/2 of the 
’’first parents”, what results in JF 2 . Then we choose a new pair of ’’parents” and repeat this 
step. The process is continued until the denominators of new ’’descendants” exceed n. Let 
p/q and p'/q' ( p/q < p'/q') be two consecutive elements of for some n. Replacing the 
pair of the ’’first parents” {0/1, 1/1} by the pair {p/q,p' / q'} and continuing the recursive 
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process described until the denominators of new ” descendants” exceed some n' > n, we 
end up with the Farey sequence of order n' generated by the ’’parents” {p/q,p' /q'}. The 
sequence so defined consists of those elements p" /q" of F n > that satisfy the inequalities 
p/q < p" j of' < p' / of . This definition can be extended to the vector case. Namely starting 
with the pair of vectors {(p/q, r), (p' /of , r')} (with p/q,p' /of as above), where 0 < r, r' < 1 
are rationals, and forming ’’descendants” componentwise we arrive at the vector Farey 
sequence of order n! (n' > n) generated by the pair of vector ” parents”. The recursion step 
is repeated until the denominators of the first components of the ’’descendants” exceed 
n! . 

For each family of the phases described so far the density pairs (p e , pf) are of the form 
(Pe/qe,Pi/q e ), where p e /q e is an irreducible rational, and the components of these density 
pairs satisfy a linear relation ap e + bp, = c. Clearly, if (p e , pf) and (p' e , p'f) are two such 
pairs, then their ’’descendant” (p'f, p'f) (defined componentwise) satisfies the same linear 
relation. Therefore the density pairs of the families of phases considered constitute the 
vector Farey sequences. Specifically the density pairs of the [p/q] n phases form the vector 
Farey sequence generated by the ’’first parents” {(0/1, 0/1), (1 /n, n/n)}, the density pairs 
of \p/q\* n phases form the vector Farey sequence generated by {(0/1, 1/1), (1/n, 0/n)}, 
while those of the \p/q]h/k and \p/q\* h /k phases form the vector Farey sequences generated 
by {(h/k, 0/k), (1/2, 1/2)} and {(h/k, k/k), (1/2, 1/2)}, respectively, provided that h/k 
and 1/2 are ’’parents” in Tk.. 

It appears that following these recursive procedures one can construct also the unit 
cells of the phases considered. For that it is enough to specify the unit cells of the 
’’first parents” and to identify the operation of forming descendants with the simple 
concatenating of these unit cells. In the case of \p/q] n phases the unit cell of the (0/1, 0/1) 
phase is an empty site with no electrons [o], while the unit cell of the (1/n, n/n) phase 

n 

consists of n consecutive occupied sites with one electron [••?.. •] . Thus, for n — 2 
their descendant has the densities (1/3, 2/3) and its unit cell is [o • •], the descendant 
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of {(0/1, 0/1), (1/3, 2/3)} has the densities (1/4, 1/2) and the unit cell [o o ••] while 
the descendant of {(1/4, 2/4), (1/3, 2/3)} has the densities (2/7, 4/7) and the unit cell 
[o o • • o • •], etc. In the case of \p/q\h/k phases the unit cell of (h/k, 0/k) phase consists 
of k unoccupied sites with h electrons and that of (1/2, 1/2) phase consists of two sites 
with one ion and one electron. 

The \p/q\h/k phases seem to be natural generalizations of the n-molecule most ho- 
mogeneous phases. They are generated just by different ’’first parents”. Therefore it is 
tempting to name the phases generated by some pair of ’’first parents” phases by means 
of the recursive concatenating the generalized most homogeneous phases. 

If on the other hand, given a pair of ’’first parents” phases, we first concatenate some 
number of the unit cells of one parent and then we concatenate the resulting composite 
cell with a similarly obtained composite cell of the other parent, we obtain a phase, which, 
compared to the corresponding most homogeneous phase, is extremely inhomogeneous. 
We term it the generalized crenel phase. The crenel phases considered in [[3j are generated 
by the unit cells {[o], [•]}. They are the ’’most inhomogeneous” phases among all the pe- 
riodic phases with given density and period. It is clear that the families of the generalized 
most homogeneous phases and the generalized crenel phases are not disjoint. 

Now we are in a position to define another family of the generalized most homogeneous 
phases, denoted \p/q\ 1 / 2 - All the members of this family have the same electron density 
p e = 1/2 and therefore they are characterized by their ion density pi = p/q. Using 
the recursive procedure of forming descendants, the density pairs of these phases can 
be obtained from the ’’first parents” {(1/2, 0/2), (1/2, 1/2)} while their unit cells are 
generated by the ’’first parents” unit cells of the form [ 00 ] with 1 electron and [«o] with 
1 electron. For instance the ion densities of the \p/q\ 1/2 phases with 1/2 < p/q < 3/4 
and whose periods do not exceed 14 constitute the following sequence (in the ascending 
order): 1/2, 8/14, 7/12, 6/10, 5/8, 9/14, 4/6, 7/10, 10/14, 3/4. The denominators of 
these fractions are equal to the periods of the corresponding phases. All these phases, 
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except those whose ion densities are 9/14, 7/10 and 10/14, are simultaneously of the 
generalized crenel type. 



The only aperiodic phases that appear in our restricted phase diagrams are those that 
constitute some mixtures of a few periodic or translationally invariant phases. The ion 
configurations of the mixtures are the corresponding mixtures of the ion configurations 
of the component phases || . A distinguished role is played by the mixtures of the trans- 
lationally invariant phases, the so called segregated phases [l|]. For p e < 1/2 only the 
mixtures of the full phases and the vacuum turn out to be stable ||. Their ion configu- 
rations are the corresponding mixtures of the full and empty configurations, i. e. all the 
ions clump together. It has been proved in ||T2[ that if p e ^ p i then there is a sufficiently 
large U such that above it the segregated phase is stable. In || the exact stability region 
of the segregated phases has been conjectured. This region consists of those points ( p e , pi ) 
that satisfy the inequality p e < Pib{U ), where b{U) is the solution to the equation derived 
in ||. More details concerning these phases can be found in [[TJ, ||, [0, |T2 . 



3 Main observations and conjectures 



The excerpts from our data are presented in Figs 1-7. First, in Figs la and lb, we show 
two schematic phase diagrams, which indicate where the families of characteristic phases, 
described in the previous section, are located in the (p e , pi)-plane. 

Then in Figs 2-6 we show the actual data in the form of the restricted canonical phase 
diagrams. The set of the admissible ion configurations used to obtain these diagrams 
consisted of all the periodic ion configurations with the period not greater than 14 and of 
the full and empty configurations. In these diagrams the (p e , pj)-phases are represented by 
the black spots. In order to determine their rational p e and pi coordinates it is sufficient to 
know that the p e values of the stable phases constitute the following sequence of irreducible 
rationals (in the ascending order): 

1/14, 1/13, 1/12, 1/11, 1/10, 1/9, 1/8, 1/7, 2/13, 1/6, 2/11, 1/5, 3/14, 2/9, 3/13, 
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1/4, 3/11, 2/7, 3/10, 4/13, 1/3, 5/14, 4/11, 3/8, 5/13, 2/5, 5/12, 3/7, 4/9, 5/11, 6/13, 
1/2. Moreover, if p e = p e /q e is any element of this sequence, then the coordinates of 
the majority of the black spots, whose p e coordinates are equal to p e /g e , belong to the 
sequence 1 /q e , 2/q e , ..., (q e — 1 )/q e - For the set of the admissible ion configurations 
used to obtain the diagrams presented the exceptions are rare (most of them are the 
phases with p e = 1/2) and they are indicated in the captions of Figs 2-6. The periodic 
ion configurations corresponding to the (p e , pi)-phases with the p e and pi coordinates 
belonging to the sequences quoted can easily be determined. Suppose that p e = p e /q e and 
Pi = p/q e . Then the unit cell has the length q e . If we label the positions of the lattice 
sites in the unit cell 0, 1, . . . , p — 1, then the positions of the ions are 0, ki, . . . , k p - 1 , 



where k\ is the inverse of p e mod q e and ki = ik\ mod q e (cf the equations [12| and the 
comments of point 2 in the sequel). 

Finally in Fig. 7 we show the intersections of the restricted canonical phase diagrams 
with the line pi = 1/2, obtained for different U values. 

1. The canonical phase diagram establishes a one to one [/-dependent correspondence 
between the density pairs (p e , p^ and the stable phases. Suppose that at the point (p e , p^, 
where p e = p e / q e and pi = Pi/qi are irreducible rationals, the stable phase is periodic. The 
stable phase has necessarily the lowest Fermi sea energy. Now according to the second 
order degenerate perturbation theory with respect to U, for given p e = p e / q e the lowest 
Fermi sea energy among the periodic ion arrangements with pi = Pi/qi is attained by 
the periodic configuration that has a gap at the Fermi level and this gap is the largest 
0- Specifically, to meet the requirement of the maximal gap, the period of the phase 



corresponding to (p e /q e , Pi/qi) has to be equal to Q - the least common multiple of q e and 



qi and the configuration w has to maximize | w (27rp e )|, where 



w 



27 r— ) = ^ exp ( i2ir—x ) w(x). 



Q - 1 



Pe 



( 11 ) 



21=0 



The same conclusion has been reached in [1]. If there is only one phase that complies 
with the rules quoted, we call it uniform (this is always the case if Q — q e ), if there is 
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more than one, we call them degenerate uniform. For all U values tested, from 0.2 up to 
6, our restricted phase diagrams conform to these rules surprisingly well. 

First, for all values of U, 0.2 < U < 6, all the phases have appeared to be insulators. 
However, in order to reach this situation for small values of U the maximal period of 
the configurations tested has to be quite large. For instance for U = 0.2, the (1/2, 1/2)- 
phase is stable also for p e in a small vicinity of 1/2, unless we take into account in the 
vicinty of (1/2, 1/2) the phases whose periods go up to about forty. Since in the canonical 
phase diagrams presented in JT|| the maximal period of the configurations tested is only 
8, for small U the phases that appear are not insulators. In the diagrams presented in 
0 this effect is even more pronounced, occurs also for U much larger than 1. This is a 
consequence of using too ’’directed” set of admissible ion configurations. 

Second, only for intermediate U values, 1 < U < 3.6 roughly, we have detected some 
nonuniform phases. One example is the (3/7, l/2)-phase whose ion configuration is not 
among those for which the maximum of | w (27rp e )| is attained. It is visible for instance 
in U = 3 phase diagram, Fig.6. Its period 14 is the least common multiple of q e and q tl 
however its unit cell is of the form [• o o • o • o o • o • • «o], while the (3/7, 1/2) phase 
that maximizes | w (27rp e )| has the unit cell [• o o • o • o • o • • o »o]. Other examples 
of nonuniform phases have been revealed by additional studies, not represented in the 
diagrams included. For U between 1 and 1.9 there appear (1/2, p^-phases whose ion 
configurations are of the generalized crenel type, generated by the ’’parents” : [o»] with 
one electron and [••] (or [oo]) with one electron. The periods of these phases are not the 
least common multipliers of q e and q t . 

2 . For every p e = p e / q e there is a characteristic U value, Ui(p e /q e ), such that below 
U\(pe/q e ) the general rule of point 1 is accompanied by an additional one. Namely, only 
the {pe/q e ,Pi/qi )- phases such that Q = q e appear in the restricted phase diagrams. The 
smallest U\(p e /q e ) is t/i(l/2) ~ 1. As a consequence, below Ui(p e /q e ) the positions of the 
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ions in the unit cells of (p e /<le,Pi/Qi ) -phases are given by the solutions kj to the equations 

p e kj = j mod q e , j = 0 , 1 ,..., p \q e / q { - 1. (12) 

It is easy to see that once the k\ solution has been found, the remaining solutions are 
given by kj = jk\ mod q e . This additional rule turns out to be satisfied also for sufficiently 
large U (above U = 3.6 roughly), independently of p e , simply because at these U values 
the ion arrangements of the periodic phases are the ion 1-molecule most homogeneous 
configurations. 

Moreover, there is another characteristic U value, t7 2 (p e /g e ), UziPe/Qe) < Ui(p e /q e ), 
such that below f/ 2 (p e /g e ), if Pi,min/qe is the minimal ion density of the stable phases with 
Pe = Pe / qe and Pi,max/qe is the maximal one, then all the remaining (p e /g e ,pj/g e ) phases 
with pi assuming all the values Pi t min, Pi, min + 1, . . . , Pi, ma x are stable as well. This is 
a more precise form of the observation made in ||. The value of U-^ipe/Qe) decreases as 
the period q e of the phases increases. It is not clear if there is some nonzero t/ 2 that 
is good for all periods. In Fig. 3 (U = 1.4), for p e = 5/14 there appear the phases 
corresponding to p t = 5, 6, 9, 10 but those with p* = 7, 8 are missing. Similarly, for 
p e = 4/11 there appear the phases corresponding to Pi = 4, 5, 7, 8 but that with Pi = 6 
is missing. Some other examples are visible in this diagram. All these missing phases 
are stable already in the phase diagram for U = 0.8, Fig. 2. The values of p ijm j n and 
Pi, max appear to be determined by the n-molecule most homogeneous phases and the 
generalized most homogeneous phases. For given p e = p e /q e the minimal Pi, m i n /q e an d 
the maximal p t . m ax/(le i° n densities are determined by the minimal and maximal pi values 
of the intersections of the line p e = p e /qe and the line segments formed by the \p/q] n , 
[p/g]n> \p/q}h/k and \p/q}* h/k stable phases. 

3. When the maximal period of the admissible ion configurations grows, then 
the number of the stable phases in the Farey sequences of the most homogeneous and 
the generalized most homogeneous phases grows as well. It appears that the order of 
any such a Farey sequence of phases is equal to the maximal period of the admissible 
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ion configurations. Moreover the Farey sequences of the most homogeneous and the 
generalized most homogeneous phases are stable in specific [/-dependent intervals of p e 
values. Only the p e intervals of \p/q) i phases are nonvoid for all U. In the case of the 
ion n-molecule most homogeneous phases with n > 1 the critical U values at which the 
p e intervals open have been determined in ||. The hole n-molecule most homogeneous 
phases appear in the phase diagrams together with the 2-molecule most homogeneous 
phases, i.e. (according to ||) below U = 2.66. The \p/q\ n families appear in the phase 
diagrams in succsession, starting with n = 2, so that for sufficiently small U only the 
phases with n not exceeding some characteristic value are stable ||. In contrast, for given 
U < 2.66 all the \p/q\* n phases with n greater than some characteristic value appear. For 
instance for U = 2 (Fig.4), those with n = 1 are still missing, but somewhere between 
[7 = 2 and U = 1.4 all the [p/g]* families become stable. The p e stability intervals of the 
\p/q]n family of phases are determined by the points of intersection of the line pi = 1 — np e 
with the segments of those lines pi = n'p e that correspond to the stable phases. 

4. The Farey sequences of the [p/q\h/k generalized most homogeneous phases and 
their adjoint phases start to appear somewhere between U — 1 and U = 0.8, as U 
decreases. In Fig. 2 (U = 0.8) three representatives of the \p/q\ 1/3 phases and also three 
representatives of the adjoint family are visible. For U = 0.4 two families, [p/q] 1/3 and 
\p/q\ 2/5 and their adjoints are well represented in our restricted phase diagrams and for 
U = 0.2 a few representatives of the \p/q] 3/7 family appear. 

5. For U between 1 and 3.6 we have found a family of phases with p e = 1/2. For 
U > 2 this family consists of the generalized most homogeneous phases of the [p/q] 1/2 
type (and the adjoint family). For instance, in Fig. 6 (U = 3) and in Fig.5 (U = 2.6) 
this family is represented by the phases whose ion densities (in the ascending order) are 
8/14, 7/12, 6/10, 5/8, 9/14 and 4/6, where the denominators are equal to the periods 
of these phases. The phases that correspond to all p. t fractions quoted, except 9/14 are 
uniform (the numerators of these pi fractions exceed the halves of the denominators by 
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one) and their unit cells are of the generalized crenel type. In contrast the (1/2,9/14)- 
phase is degenerate uniform and it is not generalized crenel. Its unit cell has the form 
[o«o»»«o»o*o* in agreement with the recursive concatenating procedure. For 
sufficiently small U the degenerate uniform \p/q] 1/2 phases that are not generalized crenel 
transform into (degenerate uniform) phases of the generalized crenel type (with the same 
ion density). Thus in some points (1/2, pi) we observe the structural phase transitions 
- one stable periodic phase transforms into another stable periodic phase as U varies. 
For instance, for U = 2.1 the (1/2, 9/14)-phase has the unit cell shown above. Al- 
ready for U = 1.9 the unit cell of the (1/2, 9/14)-phase is of the generalized crenel type: 
[o • o • o • o • o •••••] . Another manifestations of the effect described can be found in the 
diagrams for U = 2 (Fig.4) and for U = 1.4 (Fig. 3). In Fig. 4 the (1/2, pj)-phase with the 
largest ion density, pi = 7/10, is degenerate uniform and generalized crenel. Its unit cell 
has the form [o • o • o • • • ••], while the unit cell of the generalized most homogeneous 
phase [7/10] 1/2 is [o»o««»o« ••]. The remaining (1/2, pj)-phases, with equal to 
5/8, 6/10, 7/12 and 8/14, are uniform and simultaneously generalized most homogeneous 
and generalized crenel. The only (1/2, pj)-phase visible in Fig. 3, with p^ = 10/14, is also 
degenerate uniform and generalized crenel. Its unit cell is while 

the [10/14] 1 / 2 -phase has the unit cell [o • o • o • • • • • o • ••]. 

6. When U increases towards 4, all the periodic phases different from the 1-molecule 
most homogeneous ones start to dissappear (their domains in the (p e ,/q) plane shrink 
quickly). Thus, our diagrams suggest that above U > 4 the only phases that remain in 
the phase diagrams are the ion 1-molecule most homogeneous phases with arbitrary p e , 
the mixtures of the ion 1-molecule most homogeneous phases and the full phases or the 
empty phases, and the segregated phases. This is in contrast with the canonical phase 
diagrams of [|1| and 0. 
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7 . The ” evolution” of the phase diagrams with decreasing U can be viewed as 
building up an extra symmetry - the ion (or the electron) hole-particle symmetry. In 
the (p e , pj)-plane a polygon is formed, which is invariant with respect to the reflections 
in the line pt = 1/2. The segments of the lines of the ion and hole most homogeneous 
phases determined by the mutual intersections of these lines constitute the sides of that 
polygon. For instance, in the U = 1.4 restricted phase diagram (Fig.3) the tetragon with 
the following vertices (1/2, 1/2) (1/3, 2/3) (1/4, 2/4) is well visible, however it is not yet 
reflection invariant, since some phases inside this tetragon are missing, but for U = 0.8 
(Fig. 2) the reflection symmetry with respect to the line pi — 1/2 is already restored. As 
U decreases the invariant polygon gradually extends through out the smaller p e values. 
For U = 0.4 it consists of two tetragons whose boundaries are visible in the schematic 
phase diagram, Fig.l. The appearance of the hole most homogeneous phases \p/q\* n) the 
generalized most homogeneous phases \p/q\h/k and \p/q\* h / k phases, as well as the fact that, 
as U decreases the minimal p e values for which the [p/q} n phases are stable increase, can 
be viewed as manifestations of the tendency to restore the ion hole-particle symmetry. 
Because of the inversion symmetry of the phase diagram, the ion hole-particle symmetry 
is equivalent to the electron hole-particle symmetry. 

8. In agreement with the Gibbs phase rule for a two component system at a fixed 
temperature, only mixtures of two or three phases appear in the restricted canonical phase 
diagrams. Actually in the diagrams presented, Figs 2-6, there appear also mixtures of 
more than three periodic phases. This is a consequence of the finite step with which 
the electron chemical potential axis has been scanned, when the grand canonical phase 
diagrams were constructed. The stable mixtures can be arranged in five groups. For 
p e < 1/2 there appear mixtures of: (a) two periodic phases, (b) three periodic phases, (c) 
a periodic phase and a translationally invariant phase (where the electron density of the 
latter may vary in an interval), (d) two periodic phases with a translationally invariant 
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phase (whose electron density is uniquely determined by the periodic phases) and (e) 
the segregated phases that are mixtures of the vacuum (i.e. the empty phase with no 
electrons) and the full phases with p e not greater than certain [/-dependent critical value 
i, 0, 0. The mixtures belonging to the groups (a) and (b) are the insulators while 
those belonging to the other groups are the conductors. Typically the periodic phases 
that appear in the mixtures of the kind (c) and (d) belong to the Farey sequences of the 
(generalized) most homogeneous phases. When the maximal period of the admissible ion 
configurations grows, then the areas occupied by the mixtures of the kind (b) and (d) 
seem to shrink quickly. The mixtures of the kind (d), where the electrons densities of the 
two periodic phases are equal and assume a limiting value of the p e stability interval of 
one of the Farey sequences of phases, seem to be the only exceptions. In Fig. 2 these 
are for instance the mixtures of (1/6, 1/3), (1/6, 1/2) and the full phases, the mixtures 
of (2/5, 3/5), (2/5, 4/5) and the full phases, etc. Therefore we expect that in the limit of 
the infinite maximal period there are no mixtures of the kind (b) and those of the kind 
(d) are relatively rare. 

9. As U increases, then typically the periodic phases and their mixtures that 
are located at the boundary of the (p e , Pi) region, where the periodic phases and their 
mixtures are stable, lose against mixtures of other periodic phases with a full or an empty 
phase. This results in an insulator-conductor transitions. In the limit U — > oo the neutral 
phases (with arbitray p e ) remain the only insulators. However the insulator-conductor 
transitions can occur also when U decreases. As first observed in ||, as U decreases the 
n-molecular most homogeneous phases, with electron densities sufficiently close to 0 (or 
1) lose their stability against mixtures of ^/-molecular most homogeneous phases (n f > n ) 
with the empty phases. In the case n — 1 this occurs below U = 2/\/3 0. This effect is 
visible for instance in Fig.2 (U = 0.8). In the region of small p e the representatives of the 
neutral phases terminate for larger p e than this is the case for the ion n-molecular phases 
with n = 2, 3. 
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10 . If we start with the point (p e /q e ,0) and move it np along the line p e = const , 
then typically it first passes through infinitely many mixtures of some empty phases (with 
smaller p e ) and some periodic ones (with larger p e ). After that the point enters an interval 
of pi values, where only periodic phases and their mixtures are the stable phases. It is not 
clear if this interval contains only not more than q e — 1 periodic phases with the period q e 
and mixtures of two consecutive periodic phases of this kind. It may not be connected, 
there can be subintervals where mixtures of some periodic phases and full or empty phases 
are stable. Finally the point passes through infinitely many mixtures of some periodic 
phases (with smaller p e ) and some full phases (with larger p e ) until it reaches eventually 
the full phase (if p e > b(U), 0) or the segregated phases (if p e < b(U)). Similar scenario, 
but obtained by moving the point (0,1/2) along the line pi = 1/2 and then repeating 
this procedure for many values of U , is presented in Fig. 7. This kind of the restricted 
canonical phase diagram can directly be compared with the diagrams shown in |1| and 

0 . 
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Figure Captions 



Fig. 1 a) Schematic canonical phase diagram for U = 0.4. The shaded regions corre- 
spond to the segregated phases. The ion rt-molecule most homogeneous phases are 
distributed along the continuous line segments while the hole n-molecule most ho- 
mogeneous phases along the dashed line segments. The \p/q\h/k phases are lying 
along the dashed-doted line segments. 

b) Schematic canonical phase diagram for U — 4. The only periodic phases that ap- 
pear in the diagram are the ion 1-molecule most homogeneous phases. The points 
representing these phases constitute the line segment that extends from (0, 0) to 

(MJ- 



Fig. 2 The restricted canonical phase diagram for U = 0.8. This is a part of the complete 
diagram, which is sufficient for reconstructing the remaining part by means of the 
inversion symmetry (cf section 1). The black spots represent the periodic phases 
(whose period is at least 2). Those black spots whose corresponding phases coexist 
in the grand canonical phase diagram are joined by the straight line segments. The 
points lying on such a segment represent the mixtures of the two periodic phases 
that correspond to the ends of the segment. The straight line segments partition 
the phase diagram into convex areas, most of which are triangles. These triangles 
correspond to the coexistence points of three periodic phases in the grand canonical 
phase diagrams. The points located inside these triangles represent mixtures of the 
three phases that correspond to the vertices of the triangles. Some segments are 
missing and therefore there are areas which according to this restricted diagram are 
mixtures of more than three phases. This is an artefact related with the precision 
of the diagram. All the (p e , p;)-phases conform to the equations |12. 



Fig. 3 The restricted canonical phase diagram for U = 1.4. This is a part of the complete 
diagram, which is sufficient for reconstructing the remaining part by means of the 
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inversion symmetry. Only the (1/2, 10/14)-phase does not conform to the equations 

unit cell is See the caption of Fig. 2 and section 3 for more 

explanations. 

Fig. 4 The restricted canonical phase diagram for U = 2. This is a part of the complete 
diagram, which is sufficient for reconstructing the remaining part by means of the 
inversion symmetry. All the phases with p e = 1/2, except that with the largest ion 
density pi = 7/10, belong to the family of \p/q\ i/ 2 phases. Their ion densities (in the 
ascending order) are: 8/14, 7/12, 6/10, 5/8. The unit cell of the (1/2, 7/10)-phase 
is [o • o • o • • • ••]. See the caption of Fig. 2 and section 3 for more explanations. 

Fig. 5 The restricted canonical phase diagram for U = 2.6. See the caption of Fig. 2 
and section 3 for more explanations. All the phases with p e = 1/2 belong to the 
family of [p/q] 1/2 phases. Their ion densities (in the ascending order) are: 8/14, 
7/12, 6/10, 5/8, 9/14, 4/6. 

Fig. 6 The restricted canonical phase diagram for U = 3. See the caption of Fig. 2 and 
section 3 for more explanations. All the phases with p e = 1/2 belong to the family 
of \p/q\ 1/2 phases. Their ion densities (in the ascending order) are: 8/14, 7/12, 6/10, 
5/8, 9/14, 4/6. The (3/7, 7/14)-phase is not uniform. 

Fig. 7 This is a schematic phase diagram, the intersections of the restricted canonical 
phase diagrams with the line pi = 1/2, obtained for the specific U values marked on 
the horizontal axis. We can distinguish a few kinds of stable phases. In the white 
region marked ’’segregated”, the segregated phases are the only stable phases. In 
the light grey region marked ’’most homogeneous&full” a stable phase is a mixture 
of one or two ion n-molecule most homogeneous phases and a full phase. In the 
dark grey region marked ’’periodic phases and their mixtures” the stable phases are 
either periodic or mixtures of one or two periodic phases. Finally in the unmarked 



12. It does not belong also to the \p/q] 1/2 family. It is generalized crenel and its 
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grey region the stable phases can be either like in the dark or light grey regions, or 
something else, depending on p e . For the specific U values marked on the horizontal 
axis, the critical p e values, where the kind of stable phases changes, are marked 
by the short horizontal segments. The continuous segments correspond to periodic 
phases while the dotted ones to mixtures. 
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